Origin of bulk uniaxial anisotropy in zinc-blende dilute magnetic semiconductors 
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It is demonstrated that the nearest neighbor Mn pair on the GaAs (001) surface has a lower 
energy for the [110] direction comparing to the [110] case. According to the group theory and 
the Luttinger's method of invariants, this specific Mn distribution results in bulk uniaxial in-plane 
and out-of-plane anisotropies. The sign and magnitude of the corresponding anisotropy energies 
determined by a perturbation method and ab initio computations are consistent with experimental 
results. 



One of the founding blocks of condensed matter physics 
is the virtual crystal approximation [1, 2] allowing to ex- 
tend the outcome of the group theory for a given crys- 
tal symmetry to alloys with a random distribution of 
their constituents. However, it has been recently demon- 
strated, combining the progress in ab initio simulations 
and nanocharacterization methods, that the condition of 
the random distribution is violated in a number of alloys, 
leading to striking consequences. In particular, it has 
been found that open d shells of transition-metal (TM) 
cations diluted in non-magnetic compounds not only pro- 
vide localized spins but also, through hybridization with 
band states, contribute significantly to the cohesive en- 
ergy, particularly if TM atoms occupy neighboring sites 
[3|, 0] ■ The resulting attractive force between magnetic 
cations leads to their aggregation, either triggered by 
appropriate post-growth high-temperature annealing, as 
found for (Ga,Mn)As @, or occurring at the growth sur- 
face during the epitaxial process, the case of (Ga,Fe)N 
Q . The TM aggregation invalidates the main premise of 
dilute magnetic semiconductor (DMS) physics, namely 
that concerning the random distribution of TM spins over 
cation sites. One striking consequences is the appearance 
of surprising high temperature fcrromagnetism in numer- 
ous magnetically doped semiconductors and oxides, as- 
signed now to the presence of TM aggregates 0, 0| ■ 

In this paper we show that this progress in ab initio 
and nanocharacterization methods makes it possible to 
establish the origin of bulk crystalline iryjlane uniax- 
ial anisotropy found in magnetotransport [7H10( , magne- 
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12J, magnetic [12Hl5j . and ferromagnetic 
161 ] studies of (Ga,Mn)As. We show quan- 



tooptical 
resonance 

titatively that this puzzling anisotropy, whose presence 
contradicts the results of group theory for zinc-blende 
crystals, results from a non-random distribution of Mn 
over cation sites, setting in at the growth surface during 
the epitaxy. This insight allows us to propose methods 
for its controlling, the important step for further explo- 
ration of functionalities associated with its presence in 
(Ga,Mn)As and related systems |17h20| . Furthermore, 
our model elucidates the origin of a threefold enhance- 



ment of the shape magnetic anisotropy found in thin films 
of (Ga,Mn)As 

We consider here zinc-blende (Ga,Mn)As grown by 
low-temperature molecular beam epitaxy [2l[ along the 
[001] direction. Under these conditions, long-range aggre- 
gation of substitutional Mn cations is kinematically lim- 
ited, as according to three-dimensional atom probe mea- 
surements the Mn distribution is random down to at least 



1 nm [22J. Actually, the formation of Mn- rich (Mn,Ga)As 
nanocrystals inside (Ga,Mn)As films starts to be visi- 
ble under annealing at temperatures considerably greater 
than the growth temperature [HHH. Thus, we start our 
studies by finding out the energetically favorable position 
of the nearest-neighbor (NN) Mn cation dimer on the 
GaAs (001) surface, where constituent atoms are mobile 
during the epitaxy. We perform ab initio calculations 
employing SIESTA code 0, [Hj], whose localized basis is 
well suited for surface studies. The computations within 
the local-density approximation are carried out for the 
following geometry: a pair of Mn impurities is located 
on the (001) surface of a twelve monolayers thick GaAs 
supercell with lateral dimensions of two lattice constants. 
The lower Ga-terminated surface is saturated with pseu- 
doatoms of Z = 1.25, yielding a total of 112 atoms in a 
supercell, with a 1.5 nm thick vacuum region. A 5 x 5 x 1 
Monkhorst-Pack grid of fc-points is used, with the slab 
dipole correction and simulate doping options enabled. 
Since we are interested in the situation under growth 
conditions, we have performed a non-spin-polarized cal- 
culation with the lattice optimization. 

According to the computation results, the preferred 
orientation of Mn dimers is [110], with the energy gain 
comparing to the NN [110] pair being as large as 1.0 eV. 
Among many consequences of such a non-random Mn 
distribution is the appearance of a local strain, as Mn 
atoms in the pair are displaced of the GaAs cation po- 
sitions, whereas the As atom between them is shifted 
along the [001] direction by as much as 2.6% of the bond 
length. The resulting strain may contribute to the forma- 
tion of stacking faults propagating in the (111) and (111) 
planes [26| . observed in (Ga,Mn)As by high resolution 
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electron transmission microscopy (27( and synchrotron x- 
ray diffraction i28j. Importantly, according to recent ab 
initio studies [28[, the intersection lines of the stacking 
fault pairs may enhance further the aggregation of Mn 
dimers along the [110] crystallographic directions. 

We evaluate the effect of this specific distribution of 
Mn ions upon magnetic anisotropy in two steps. First, by 
employing the group theory and the Luttinger's method 
of invariants we establish the expected form of magnetic 
energy depending on the magnetization orientation. Sec- 
ond, by making use of either perturbation theory or an ab 
initio method we evaluate the magnitude of anisotropy 
constants brought about by the non-random Mn distri- 
bution. 

More specifically, we place a nearest neighbor (NN) Mn 
pair in a GaAs supercell. In order to consider various 
spatial correlations between these pairs, the supercells 
are arranged in two ways, corresponding to positions of 
the Mn dimers in a simple cubic or face-centered cubic 
sublattice, respectively. Such a system possesses the C 2 „ 
symmetry. If the Mn pair is along [110] direction, this 
group comprises a twofold rotation axis [001] and a reflec- 
tion plane (110). Accordingly, the magnetic anisotropy 
energy F(m x , my, m z ) as a function of the magnetization 
direction (rn x ,m y ,Tn z ) — (sin cosy, sin $ sin ip, cos &) 
can be expanded into invariants fi, 



A 



(1) 



where the basis invariants are defined respecting the de- 
composition of the space of spherical harmonics with 
given I into irreducible representations of the Td group. 
For example, the decomposition E © T2 of the space 
of spherical harmonics with I = 2 defines invariants 



results in a non-zero value of the macroscopic parameter 
A'2. The macroscopic value of A2 is a product of the 
anisotropy energy of a single Mn pair and the difference 
in the densities (per supercell volume) of pairs along the 
[110] and [110] directions, respectively. 

More formally, the correlations in occupations between 
the neighboring sites form a representation of the point 
group Td j which decomposes into a sum of the irreducible 
components, A\ © E © T2. The irreducible components 
are the total pair density (A± ) , the distribution of pairs 
among the three planes (E) 1 and the differences be- 
tween occupations of the two directions in each of the 
planes (I2). Similarly, the magnetic anisotropy free en- 
ergy function decomposes into E © T2 (for I = 2) and 
Ai © E © Ti © T 2 (for 1=4). Here, we consider a first 
order (linear) dependence of the magnetic anisotropy on 
the pair correlations. Hence, the irreducible components 
are in a direct correspondence: the cubic anisotropy pa- 
rameter depends linearly on the total pair density (rep- 
resentation A\), the out-of-plane uniaxial anisotropy pa- 
rameter K\ on the distribution of pairs among the planes 
(A), and the in-plane uniaxial anisotropy constant A2 
on the difference between the densities of the [110] and 
[110] pairs (T2). The remaining component, Ti, has no 
counterpart in the decomposition of the correlation func- 
tion and therefore has to vanish. We conclude that the 
non-random distribution of Mn cations of the form in- 
troduced here leads to the functional of the magnetic 
anisotropy e nergy (E q. [1]) consistent with experimental 

results dEMIll!]- 

Previously, the anisotropy constant A2 accounting for 
experimentally observed bulk uniaxial in-plane crystal- 
lographic anisotropy was described theoretically by as- 
suming the presence of shear strain, whose magnitude 



fi 



i (representation A) and / 2 = m x m y (rep- e xy was treated as an adjustable parameter [lj, |29|, |3C 



resentation T2). These are out-of-plane and in-plane 
uniaxial anisotropies. For 1 = 4, the decomposition is 
A\ © A © Ti © T2, where the invariant from A\ is the cu- 
bic anisotropy, the ones from A and T2 are higher order 
uniaxial anisotropies analogous to those with I = 2, and 
there is no invariant in T±. 

Let us consider now possible NN cation configurations. 
Up to the translational symmetry, there are six different 
placements corresponding to the three planes [i.e., (100), 
(010), and (001)] in which a pair can be located and two 
possible directions in each of the planes. As the growth 
occurs along the [001] direction, we cannot expect the 
same density of dimers in the (001) plane comparing to 
the (100) and (010) planes. Thus, averaging of f± over 
the three planes may lead to K\ ^ 0, adding to the ef- 
fect of epitaxial biaxial strain. Furthermore, by symme- 
try, the densities of NN pairs in the two directions in the 
(100) and (010) planes are equal, therefore / 2 averages to 
zero for each of those two planes. However, the symme- 
try allows a different density of NN pairs along the [110] 
and [110] directions in the (001) growth surface. This 



evaluated with the p — d Zcner model to be of the order 

of e 
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0.05% for the relevant GaAs deformation po- 
tential d = —4.8 eV 2iJ 30| . Such a shear deformation 
of (Ga,Mn)As films has not been found experimentally 
[28l | . Similarly, an additional contribution to K\, found 
in studies of magnetic anisotropy as a function of biaxial 
strain [101] . can be parametrized by e xx « —0.05%. We 
will show now that new terms in the Luttinger hamilto- 
nian brought about by the non-random Mn distribution 
have the form of a strain hamiltonian whose elements 
have signs and magnitudes consistent with the experi- 
mental values. 

We first examine how lowering of symmetry to Civ af- 
fects the three-band effective mass hamiltonian describ- 
ing the valence band in the cubic case. The spin-orbit 
and p — d interactions can then be taken into account 
in the standard way leading to the six-band hamiltonian 
from which the magnetic anisotropy energy A (Eq. [1]) 
can be directly determined [29L l30f. By employing the 
Luttinger's method of invariants [311 ] . we find that the 
the presence of spatial correlations of Mn ions is cap- 
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FIG. 1. (Color online) Effective shear e xy (a) and epitaxial 
biaxial e xx (b) strains, normalized to unit pair density, for 
a nearest-cation neighbor [110] pair of Mn ions in supercells 
arranged in simple cubic and face-centered lattices, as calcu- 
lated by using Eq. (j4|). 



tured by the use of the virtual crystal k ■ p hamiltonian 
with terms corresponding to effective shear and biax- 
ial strains, described by two components of the strain 
tensor, e xy and e xx , respectively, together with the cor- 
responding deformation potentials d = —4.8 eV and 
b = —2.0 eV. The effective biaxial strain e xx (with 
£ yy = s xx , £zz = -2ci 2 s xx / Cn, C\ijc\\ = 0.453) will 
renormalize the magnitude of strain coming from a mis- 
match to the substrate. 

In order to evaluate the magnitude of e xy and e xx , we 
start from the three-band unperturbed hamiltonian, 

h 2 

^3x3(k) = E v l 3x3 



H, 



X 



2m 

Akl+B{k 2 y +k 2 z ) Ck x k v 

Ck x k y Ak 2 y +B(kl+kl) 



(2) 



Ckxk z 
Akl+B{kl+k 2 v ) 



where in terms of the Luttinger parameters 7^, A = 
71 + 472, B — 71 — 272, and C = 673. The Mn poten- 



tial is assumed to contain a contribution from a screened 
Coulomb part and a central cell correction, leading in the 
fc-space to 



V(k) = - 



1 



ee k 2 



5/4 



(3) 



where the static dielectric constant e = 12.9, the screen- 
ing radius is r s = 0.5 nm, and the Gaussian central cell 
correction parameters are r$ = 0.28 nm, Vq = 3.0 eV [32j |. 

Within the second order perturbation theory for a Mn 
pair we obtain at k = 0, 

2 



AH. 



(2) 
3x3 



3d 



d 3 k 



(2tt) 3 

X [JS„l3x3 



2 cos ■ 



kd 



\V(k)\ 



#3x 3 (k)] 



(4) 



where d is the vector between the two Mn ions forming 
the pair and is the density of such pairs. In the case 
of a supercell, the integral should be replaced with a sum 
over the reciprocal lattice J (f^p ~ ^ y- Ek^c where 
k = is omitted from the sum and V sc is the volume of 
the supercell. The above formula assumes additivity of 
the Mn potential and neglects distortion of the lattice in 
the presence of the Mn pair. A similar formula without 
the squared cosine factor can be used to determine the 
anisotropy of a single Mn acceptor in a non-cubic super- 
cell [HI], as a single Mn acceptor residing in a supercell 
exhibits a significant magnetic anisotropy in accord with 
the symmetry of the supercell. 

The obtained results are presented in Fig.[T]for the NN 
Mn pair residing along the preferred direction [110]. In 
particular, in the case of no supercell, that is when posi- 
tions of the [1 10] dimers can be regarded as uncorrelated, 
we obtain £ xy /gd — 0.019 nm 3 , e xx /g<i — — 0.0015 nm 3 . 
For the Mn concentration x = 6.25%, this leads to 
e xy — 1.33%, e xx = —0.10%. However, when the po- 
sitions of the [110] dimers are correlated, e xy = 0.77%, 
e xx = —0.16% in the case of a 54 atoms fee supercell and 
£ X y = 0.70%, e xx = —0.26% in the case of a 64 atoms sc 
supercell, corresponding to (Ga,Mn)As with 7.41% and 
6.25% of Mn arranged into a regular sublattice of [110] 
Mn dimers. 

Comparing to experimental findings, we see that the 
computed values of e xy and e xx have the correct sign. 
At the same time, their absolute values are significantly 
greater that the one determined experimentally, e xy ss 
0.05% and e xx s» —0.05%. This could be expected as in 
real samples only a fraction of the Mn content forms NN 
pairs as well as some Mn pairs choose the less preferred 
direction during the growth process. 

It is instructive to compare the above findings to ab 
initio result. For this calculations particularly suitable 
is the Quantum ESPRESSO code developed in the plane 
wave basis 34[ . Within local spin density approximation 
and for supercells with 54 cations (fee) and 64 cations 
(sc), using a 4 x 4 x 4 grid of fc-points, we obtain from 
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magnitudes of the valence band splitting at the T point of 
the Brillouin zone, e xy = 1.70% (1.79%), e xx = -0.75% 
(-1.01%) and e xy = 1.37% (1.50%), e xx = -0.49% 
(—0.47%) without (with) lattice optimization, respec- 
tively. 

In summary, we argue that puzzling bulk in-plane 
crystalline magnetic anisotropy in (001) (Ga,Mn)As is 
brought about by the preferred formation of Mn dimers 
along the [110] at the growth surface, as implied by the 
ab initio results. The group theory and the Luttinger's 
method of invariants as well as the perturbation and ab 
initio computation results show that the effect of the pre- 
dicted Mn distribution can be parametrized by effective 
shear and biaxial strains. Their signs agree with the ex- 
perimental determination while the computed absolute 
values are much larger than the experimental magnitudes 
indicating that the surplus of Mn dimers residing in the 
preferred positions is only partial after the growth pro- 
cess. This suggests that it might be possible to control 
the strength of uniaxial anisotropy by changing the epi- 
taxy conditions, particularly the growth rate and/or tem- 
perature. Furthermore, when reducing the film thickness, 
intcrfacial and surface anisotropies, for which in-plane 
uniaxial anisotropy is allowed even for a random distri- 
bution of magnetic ions 35 , 3(| , may gradually come into 
play. In general terms, our results show how a specific 
microscopic distribution of alloy constituents may affect 
the symmetry properties and magnitudes of macroscopic 
response functions. 
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